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ABSTRACT 

Recent radio observations reveal the existence of mini radio lobes in active galaxies 
with their scales of ^ 10 pc. The lobes are expected to be filled with shock accelerated 
electrons and protons. In this work, we examine the photon spectra from the mini 
lobes, properly taking the hadronic processes into account. Wc find that the resultant 
broadband spectra contain the two distinct hadronic bumps in 7-ray bands, i.e., the 
proton synchrotron bump at ^ MeV and the synchrotron bump at ~ GeV due to the 
secondary electrons/positrons produced via photo- pion cascade. Especially when the 
duration of particle injection is shorter than the lobe age, radio-dark 7-ray lobes are 
predicted. The existence of the 7-ray lobes could be testable with the future TeV-7 
telescope CTA. 
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1 INTRODUCTION 

Thanks to the progress of VLBI (Very Long Baseline Inter- 
ferometry) observations, compact radio lobes with a linear 
size LS ^ 1 kpc, defined as a projected length from the core 
to the lobe, have been discovered (e.g., Fanti et al. 1995; 
Readhead et al. 1996; O'Dea & Baum 1997). Further VLBI 
observations recently reveal the existence of very small radio 
lobes with LS 10 pc among two samples of compact ra- 
dio sources. The sample of compact radio sources with their 
spectral peaks higher than 5 GHz is termed as high fre- 
quency peakers (HFPs). Some of them turn out to be mini 
lobes based on their morphologies and non- variabilities (e.g., 
Orienti et al. 2007; Orienti & Dallacasa 2008) . The sample of 
compact radio sources at low-redshift {z < 0.16) is selected 
and called as CORALZ. The VLBI observations of CORALZ 
show that some of them are also found as mini radio lobes 
(Snellen et al. 2004; de Vries et al. 2009). Apart from the 
above two cases, recurrent radio sources are also known to 
possess mini lobes inside large lobes (e.g., 3C 84, Walker et 
al. 2000). These mini radio lobes are young and their ages 
are typically estimated as tago ^ 10^""^ yrs (e.g., Fanti 2009, 
Giroletti & Polatidis 2009). High-energy emission of the mini 
radio lobes have been theoretically explored by some authors 
(Stawarz et al. 2008; Kino et al. 2007, 2009). However, pre- 
vious work has focused on leptonic processes and little is 
known about hadronic processes in the lobes. 

In this Letter, we indicate that the hot spots in the 
mini radio lobe can be a plausible site of proton acceleration. 
Therefore, high energy protons are naturally expected in the 



mini lobes. In order to constrain the amount of high energy 
protons, next we examine the predicted photon spectra. Be- 
cause of their smallness, observed mini radio lobes generally 
have dense synchrotron photon fields. Furthermore, the mini 
lobes are close to the AGN core and illuminated by the core 
emission. Therefore p7 interaction is inevitable for the mini 
lobes. We calculate the photon spectra from the mini lobes 
taking into account the hadronic processes and show that 
the spectra are useful for constraining the amount of pro- 
tons via 7-ray emission features. The predicted spectra will 
be useful for testing whether cosmic ray acceleration indeed 
takes place in mini lobes. 



2 HADRONIC EMISSION 

When a jet interacts with the surrounding medium, most of 
its kinetic energy is dissipated via shocks. The termination 
point of the jet is known as a hot spot and it is identified as 
the reverse shocked region. The shocked jet plasma escapes 
sideways from the spots forming a pair of radio lobes. There- 
fore, radio lobes are remnants of the decelerated jets which 
contain relativistic particles (e.g., Begelman et al. 1984). As 
for a ratio of electrons to protons in AGN jets, it is still 
under debate (e.g., Kino and Takahara (2004), Birzan et al. 
(2008)). There are case studies indicating that protons are 
dynamically dominant at hot spots of radio galaxies (e.g., 
Stawarz et al. 2007). Here we examine what happens when 
assuming the lobe contains relativistic protons. The notation 
Q = Qx X 10^ is used unless otherwise noted. 
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2.1 Proton acceleration at the hot spots 

Hot spots in powerful radio galaxies are one of the most 
promising sites for proton acceleration (e.g., Rachen and 
Biermann 1993). Here we show that hot spots in mini lobes 
are also possible sights for proton acceleration. The energy 
of shock accelerated protons is given by Sp = "fpnipC^ where 
7p is the proton Lorentz factor measured in the shock frame. 
Although a Bohm-type diffusion timescale is a rough approx- 
imation, we adopt it for simplicity. Then the acceleration 
timescale at the hot spot fp,acc,hs = ip£p/ (eShsC) in the case 
of relativistic shocks can be estimated as 



3.5 £p,i8Bhs!-iCp.2 yr. 



(1) 



where Bhs, and are the magnetic field strength, and the 
Gyro-factor for proton accelerations at the hot spot, respec- 
tively. Although there is no direct estimate of Bhs, magnetic 
field strength in mini lobes (Biobo) has recently been esti- 
mated as Biobc ~ 10 — 100 mG by VLBI observations (Ori- 
enti & Dallacasa 2008). Since the magnetic field strength at 
the hot spot Bhs should be comparable to or larger than 
Biobc, here we set Bha ~ 100 mG. The value of ^p is a free 
parameter. In the case of blazars, ^p is not expected to be 
smaller than 10 (Inoue and Takahara 1996). As for mini 
lobes, there is little constraint on ^p. As a first step, we only 
treat the case of ^p = lO'^ in this Letter. The case of ^p < 10^ 
is not observationally excluded and it is worth to examine it 
in the context of highest energy cosmic ray sources. When 
3> 10^, photo-pion cascade does not occur. So, we do not 
treat it. The timescale of proton synchrotron cooling is 



tr. 



4„3 



,syn,hs 



2, R2 1-4 X lO"* epjgBhB^i yr, 



(2) 



where ax is the Thomson cross section. The escape velocity 
of shocked matter from the hot spot via sideways expansions 
is typically « 0.3 c (e.g.. Kino and Takahara 2004). 

Then the escape timescale is foschs ~ 3 i?hs,i8 years where 
Rhs is the hot spot radius. Since fcsc.hs "C fsyn.hs is satisfied, 
the maximum energy of protons Ep.max is obtained by the 
relation fp.acc.hs = ^csc.hs and it is 



0.9 X 10'*Bhs,-ii?hs,i8ep:2 eV. 



(3) 



The Larmor radius of protons in the spot (^L.hs) satisfies 
Rhs > ri^,hs ~ 1 X 10~^ EpAS^ha-i PC- Thus we find that the 
hot spots in mini lobes are feasible sites for proton acceler- 
ation. The escaped protons from the spots are then injected 
into the mini lobes and subsequently undergo cooling there. 



2.2 Proton cooling in the mini radio lobes 

The relativistic protons injected in mini lobes undergo 
various coolings. Before we describe the results of the 
Monte Carlo simulation in detail, it is worth estimat- 
ing the relevant cooling timescales. First, the adiabatic 
loss timescale at the mini lobes is given by tad ~ 
65 (i?iobe/2 pc)(wcxp,iobc/0.1c)~^ yrs where Riahe and 
fexp,iobe are the radius and expansion velocity of the lobe, re- 
spectively. As shown above, the Larmor radius of protons of 
10^* eV is sufficiently smaller than -Riobe- Therefore, they 
are confined in the lobes and suffer adiabatic expansion 
loss. Next, we estimate the timescale of the photo-meson 



loss process of high energy protons. As shown in the in- 
troduction, p7 interactions are expected in the lobe where 
target photons are UV photons from the accretion disk and 
synchrotron photons emitted by the lobe itself. The con- 
dition to create pious is epe-,,i ^ 0.2 GeV^ where El,i IS 
the target photon energy (i — disk, syn). Then, the typical 
energy of protons interacting with the accretion disk pho- 
tons is given by £p « 2 x 10^®(e^,disk/10 eV)"^ eV. The 
timescale of thep7 interaction is tpj « {nj^iapjC)~^ ~ 1.2 x 
lO" (LS/10 pc)~^LdiBk, 45(^7. disk/lO eV)~^ yr where cTp^ ^ 
5 X lO"'^* cm~^ at the resonance peak (Waxman and Bahcall 
1997) and Ldisk and n-y^disk are the disk luminosity and the 
photon number density n-y,disk ~ I/disk/(47rLS'^ce-y,disk), re- 
spectively. For I/disk, we assume a value of observed big blue 
bumps in Seyfert 2 galaxies (e.g., Koratkar and Blaes 1999). 
In this case, we have tad/tp7 ~ 5.3 x 10^''-Riobci<S'^i(^isk ^ 
Riahc when we assume LS oc -Riobc. For protons with Sp > 
2 X 10^® eV, the target photons are the electron synchrotron 



,/(^i?i 



■lobe'-<^7>syn^ 



ones. In this range, n.^,syn ~ 

1.7 X 10-* Lsyn,44(i?lobc/2 pc) (£^,eyn/10-'' eV)"' Cm'^* 

where Lsyn is the synchrotron luminosity emitted by 
primary electrons. Then we have tad/tp7 ~ 5.3 x 

10-* i?fobe4yn,44(e7.Byn/10-« cV)"^ CX <bc- So, tp-, bc- 

come shorter than tad when 7?iobc is sufficiently small, al- 
though such a case is beyond the scope of this work. The 
timescale of proton synchrotron loss at the lobe has been 
already given by Eq. since we assume Biobo ~ Bhs for 
simplicity. 

The timescale of energy loss by proton-proton collisions 
is given by tpp ~ 5 x 10'^(next/l cm"'')-^ yr (Sikora et al. 
1987) where ricxt is the number density of the ambient mat- 



ter. Here we focus on the case of tpp 2> ta 



10^ 



' yr. How- 



ever, there is a possibility that mini lobes expand in rather 
dense environments involving clouds of gas emitting narrow 
emission lines with next 3> 1 cm^'', and that tpp ~ tagc 
holds. We will examine it in our future work. 



3 RESULTANT SPECTRA 

To calculate broadband photon spectra, the Monte Carlo 
simulation has been performed in this work. The details of 
numerical code have been described in Asano et al. (2007, 
2008, 2009) and references therein. Here we briefly review 
the code. We calculate steady state spectra of protons {i = 
p) and electrons (i — e) by solving the following kinetic 
equations 



qi{ei) = --^ [ni{e^)ei] + — , 

OEi tad 



(4) 



where ni{ei), ii, and qi{ei) = Kte^ "' exp ^ ^ . ^^^^ ^ for 

Ei.min ^ £i are the number density, the cooling rate, and the 
injection rate of relativistic particles, respectively. Distribu- 
tions of particles and photons are assumed to be isotropic. 
We assume that the common Fermi acceleration process 
works both for electrons and protons. In this case, the re- 
sponse to magnetic turbulences is same for relativistic par- 
ticles with a same energy. Therefore, it is conservative to 
set Se = Sp and = Cp- As explained in 2.1, ep,max are 
uniquely determined for given hot spot parameters. The 
maximum energy of electrons can be obtained as ee,max ~ 
10 B^yjj^£_~l^'^ TeV from the relation te,syn,hs = te,acc,hs 
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where te,ayn,hs and te ,acc,hs " (Ce£e)/eBhsC are the syn- 
chrotron coohng and acceleration timescales for electrons, 
respectively. The kinetic equation of photons with the en- 
ergy is given by 



(5) 



where A'^-y(e-y) is the total production rate of each population 
of photons and n^(e^) is the corresponding photon number 
density. As for A''^, and 7p, we include the following physical 
processes: (1) photo-pion production from protons and neu- 
trons {p + j — >■ p/n-t-TT^/TT^), (2) the pions decay (tt" — ;> 27) 
and [tt'^ -> /i^ + i^^j —5- e+ + !/e + i^^j + (3) photon-photon 
pair production (7 + 7 — !> e"'" + e^), (4) Bethe-Heitler pair 
production (p + 7 — !> p + e"'" + e^), (5) synchrotron and in- 
verse Compton processes from electrons/positrons, protons, 
pions, muons with Klein-Nishina cross section, and (6) syn- 
chrotron self-absorption for electrons/positrons. 

The model parameter values are summarized here. For 
the hot spot, we set Rhs = 0.3 pc, Bhs = 0.1 G, and 
'^csc.hs ~ c/3. These hot spot quantities uniquely determine 



We adopt — 



1 X 10^ 



2, and 



ep,min/"ipC = ee,min/"ieC = 10 for particlc injections. The 
mini lobes parameters are -Riobe = 2 pc, -Biobe = 0.1 G, 
fexp,iobe = 0.1 c, aud LS — 10 pc. The accretion disk quanti- 
ties are assumed as Ldisk = 3x 10^^ erg s^^ with the average 
temperature 10 eV. Because of the obscuration by surround- 
ing dusty torus Ldisk in mini radio lobes is not directly con- 
strained. Following the work of Ostorero et al. (2010), here 
we assume Ldisk as the same one in luminous blazars. The in- 
jection power of protons Lp should be smaller than the total 
kinetic powers of powerful jets ~ ]^q47-48 g-i j^^q et 

al. 2008, GhiseUini et al. 2009). Here Lp = 5 x 10*® erg s"^ 
is assumed. Various injection powers of electrons Le will be 
examined below. 



3.1 Leptonic-hadronic model 

Fig. 1 displays the resultant photon spectrum for = 1 x 
10*^ erg s~^ (the thick solid line). There are several distinct 
features in the spectrum. 

(i) The bump at ~ sub PeV energies is composed of tt"- 
decay photons. The 77 absorption opacity against the target 
photons with the energy e7,tgt at the threshold £767, tgt ~ 
(2meC^)^ is written as T^j[e~/) = (3/8)crTn-y,tgt-Riobc (e.g., 
Razzaque et al. 2004) and it is estimated as 

Rlobe 



y{e^ = 10''' eV) 



0.8 



3 X 10« erg s-i 2 pc 

c- \ -1 

t-y,tgt 



. 10-2 eV 



(6) 



where the number density of target photons is n^^tgt = 
N-y,tgt/{nRf^i,^c) = L7,tgt/(7re7,tgt-R?obeC). Hence, the 7r°- 
decay photons with e-y ~ 10^* eV can partially escape from 
the lobe. The escaped 100 TeV photons, however, would 
be absorbed by the radio background photons during their 
propagations since their mean free path is shorter than 
10 Mpc (e.g., Coppi and Aharonian 1997). 

(ii) The secondary electrons/positrons produce the syn- 
chrotron bump at ~ GeV. They are produced via fi- and 
TT-decays. The pairs created by the 77 absorption of the 



tt" -h- 27 photons have an energy of 1/2 x 1/2 = 1/4 that of 
the parent tt". In the tt^ decay mode, the comparable frac- 
tion of the energy is converted from the parent vr"*" to the 
positron. Then, the synchrotron emission from secondary 
electrons/positrons with ~ 10^^ eV has its peak at ~ GeV. 

(iii) The proton synchrotron bump (p-SYN) appears 
at MeV energy band with the peak at i'p,syn ~ 
(epeBiobo/27rm^c^) ^ 0.7 ep.igSiobc.-i MeV. Below ~ keV, 
the synchrotron emission from the primary electrons over- 
whelms the proton synchrotron emission. Protons also trans- 
fer energy to photons by inverse Compton scattering and the 
photon spectrum is shown in the thin line (p-IC). However, 
they do not escape from the lobe because of the 77 absorp- 
tion. 

(iv) The break of synchrotron spectrum at ~ GHz is 
well known synchrotron self absorption (SSA) observed 
by VLBI observations (e.g., Snellen et al. 2000). The 
SSA turnover frequency for Se = 2 is given by Vssa ~ 



0.39 B, 



1/7 



2/7 



-4/7 



2 pc 



GHz (Keller- 

mann and Pauliny-Toth 1981) where L^yn > Lssa = 
I'ssa-Li^saa ■ The electrons with 7e ~ 100 emitting synchrotron 
at ~ 1 GHz will cool down, since the synchrotron cool- 
ing timescale satisfies te,syn ~ 25 Biobo,-i(7e/10')Vr < 
tad. Therefore the cooling break frequency obtained from 
ie.syn = tad is Comparable to or below Vssa- This implies 

that Le ~ Lsyn • 

(v) For comparison, the case of no-proton injection (i.e., 
Lp = 0) is shown in the thin solid line. By comparing the 
spectrum to the one with proton injection, we can clearly 
recognize the contribution of the hadronic emission in the 
7-ray energy domain. 

Fig. 2 shows the photon spectra for Le = 1 x 10*^ erg s"^. 



1 X 10"" erg s-\ 
1 X lO"*' erg s-^ 



1 X 10*^ erg s"\ 



1 X 10*^ erg s-\ 



and 

The Lp/ Le ratio in AGN jets is not well 
known. The examined Lp/Le ratio here partly exceeds the 
ratio of ~ 10^ measured around the Earth. However, when 
a jet with Lp/Le ^ 10^ gradually decreases its power keep- 
ing the Lp/Le ratio, the resultant spectrum may be similar 
to the one for larger Lp/Le because of fast cooling of elec- 
trons. The spectra with and without proton injection are 
displayed with thick and thin solid lines, respectively. For 
Le = 1 X 10*^ erg s~^, the leptonic emissions overwhelm 
hadronic ones at all energy domains. As Le decreases to 
Le ^ 1 X 10'*'* erg s"^, thick and thin lines become separa- 
ble in the 7-ray domain because the synchrotron emission 
from the secondary electrons/positrons overwhelms the in- 
verse Compton component from primary accelerated elec- 
trons. In MeV range, the proton synchrotron bump appears 
for smaller Le. Typical Lgyn depends on the source pop- 
ulation. CORALZs typically have Lssa ~ iO'*f-42 gj.g g-i 
while HFPs have Lssa ~ ig**"*® erg s^'. Observed turnover 
frequencies of mini lobes shows i^ssa ~ 0.1 — 10 GHz. The 
maximum ^ssa in Fig. 2 is smaller than the ones for HFPs, 
mainly because the angular size of the lobe in our model is 
larger by a factor of a few than the one for HFPs. It does 
not affect the main results of this work. The synchrotron self 
Compton component appeared in 7-ray bands is less dom- 
inant than the synchrotron one since Usyn ^ 



4 X 10 -Bfobe -1 erg cm holds for Ls; 



< lO*'^ erg s-^ 



where f/syn is the energy density of synchrotron photons. 
When Le < 10** erg s"*, the TeV photons in the lobe be- 



4 M. Kino, K. Asano 



come transparent against 77 absorption because of r^.^ (e^ = 
10^2 eV) ~ 0.8(L^,tgt/3 x lO'*^ erg s-^)(e^,tgt/l eV). Then 
the TeV photons begin to escape from the lobe. 

3.2 Pure hadronic model 

In Fig. 3, we show the spectrum with the same parameters 
as Fig. 1 but for pure proton injection (i.e., Le = 0) with 
the proton injection duration finj = 90 yrs. In order to re- 
alize Le = approximation, tinj < tagc and te,ayn < iagc 

should be satisfied where fe.syn is the synchrotron cooling 
timescale for electrons. In this case, the injection has been 
already stopped and primary electrons have been already 
cooled down. Since the source age in this model can be es- 
timated as fago ~ LS/vexp ~ 3 X 10^ yr and it agrees with 
the observational estimates of tagc ~ 10^""^ yr (e.g., Giro- 
letti 2009), the conditions of tinj < tagc and fage > te,ayn ~ 
25 -BiobG,~i(7e/10^)^ yr are indeed justified. The model con- 
sidered here is not applicable when B satisfy tage < te,syn- 
Short tinj would be naturally realized for jets with inter- 
mittent activities. For example, the recurrent mini lobe 3C 
84 actually show tinj ~ 50 yr (Asada et al. 2006). The low 
power mini lobes which seem to be dying ones (Giroletti et 
al. 2005) could also be candidates. In this model, the de- 
posited energy of protons is given by Lptinj. 

The prominent three bumps of hadronic emissions, i.e., 
the 7r°-decay photon bump at ~ PeV, the synchrotron one 
from secondary electrons/positrons at ~ GeV, and the pro- 
ton synchrotron one at ~ MeV, which have also been pre- 
dicted in Fig. 1 emerge more clearly. The escaped PeV pho- 
tons from the lobe are absorbed by the radio background 
photons and they may not reach the Earth. 

Below the optical energy band, the synchrotron emis- 
sion from secondary pairs via Bethe-Heitler process 
(P7 —7- pe~e'^) is responsible for the emission. The protons 
with Ep ~ 10^* eV interact with the disk photons and con- 
vert a fraction of ~ nie/mp of their energy to create 
pairs. Then the pairs with ~ 10^^ eV radiate synchrotron 
emission peaked at ~ optical band. Since the lobes are dim 
in radio band, we may call them radio-dark 7-ray lobes. 



cooled down but protons have not. Then, the predicted emis- 
sion is purely hadronic and the hadronic bumps are clearly 
seen. Importantly, the high energy tail of the GeV bump 
is detectable by The Cherenkov Telescope Array (CTA) 
(http://www.cta-observatory.org/,). These sources may be 
identified as radio dark 7-ray lobes. 

(iii) The predicted X-ray flux is well above the detection 
limit of XMM. The observations actually show bright X- 
ray emission. This is traditionally interpreted as thermal 
radiation from the accretion disk and the possibility of lobe 
emission has been alternatively indicated (Ostorero et al. 
2010 and reference therein). In any case. X-ray emission is 
likely composed of various different components. Therefore, 
it seems difficult to extract the hadronic component from 
the X-ray band. 

(iv) We comment on the importance of larger next it may 
lead to frre-free absorption (FFA) in radio band (e.g., Begel- 
man 1999; BickneU 2003; Stawarz et al. 2008). Actually, the 
low frequency turnover in the radio spectra of some sources 
are indeed reproduced by FFA and not by SSA (e.g., OQ 
208, Kameno et al. 2000; 0108-1-388, Marr et al. 2001). It is 
clear that dense emvironments lead to an effective proton- 
proton collision. We will examine it in the future. 

(v) We add a comment on the recent VLBI observation 
of mini lobe 3C 84. It shows the outburst around 2005 
and a new component smaller than 1 pc emerges (Nagai 
et al. 2010). Fermi/LAT also detect GeV 7-ray emission 
from it (Abdo et al. 2009). A future collaboration with 
Space VLBI project VSOP-2 with high angular resolution 
( http://www.vsop.isas.jaxa.jp/vsop2/ j) and CTA would pro- 
vide us valuable constraints on the hadronic model. Theo- 
retically, we plan to conduct studies with smaller Riohe in 
our future work. 
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4 SUMMARY AND DISCUSSION 

In this work, we point out that the hot spots in mini lobes 
are feasible sites of proton acceleration. Next, the expected 
photon spectra of the mini lobes including the hadronic pro- 
cesses are explored. Summary and discussions are as follows. 

(i) For bright lobes with Le ~ isyn ~ 10** erg s~^, the 
predicted high energy emission is detectable with the cur- 
rent 7-ray telescopes. However, it is overwhelmed by the 
leptonic inverse Compton component, and it seems hard to 
test whether the emission is of hadronic- or leptonic-origin. 
For Le ~ Lsyn ^ 10'*'' erg s ^, proton synchrotron bump 
appears at ~ MeV and the synchrotron emission from the 
secondary electrons/positrons generated by the /i- and tt- 
decays emerges at ~ GeV. The two distinctive bumps in 
7-ray domain are of hadronic origin. 

(ii) The case of the short term particle injection is ex- 
amined. It may be realized when jets have intermittent ac- 
tivities. Typically, primary accelerated electrons have been 
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Figure 1. Predicted photon spectrum from tlie mini radio lobe 
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3 X 10^^ erg s~^. The thick solid line represents the spectrum 
including the effect of internal 77 absorption. Thin blue lines 
show the proton emissions, while thin red lines represent electron 
emissions before including the 77 absorption. The luminosity is 
converted to the flux with the distance 100 Mpc. 
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Figure 2. Same as Fig. 1 but for various Lc from Lc = 1 x 
10*5 erg s~^ to 1 X 10*^ erg s~^. The sensitivities of Fermi/LAT, 
and HESS are, respectively, the ones for 1 year, and 50 h integra- 
tion time. 
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Figure 3. Same as Fig. 1 but for the case pure proton injection 
*inj = 90 years. The synchrotron emission from electrons/protons 
produced via Bethe-Heitler process and pion- and muon-decay are 
predicted. The Proton synchrotron peaks at ~ MeV. PeV pho- 
tons via pion-decay are absorbed during their propagations. The 
sensitivity for XMM-Newton is normalized by 100 ks integration 
time. 



